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Abstract. Cluster media are dynamical, not static; observational evidence suggests they 
are turbulent. High-resolution simulations of the intracluster media (ICMs) and of ideal- 
ized, similar media help us understand the complex physics and astrophysics involved. We 
present a brief overview of the physics behind ICM turbulence and outline the processes 
that control its development. High-resolution, compressible, isothermal MHD simulations 
are used to illustrate important dynamical properties of turbulence that develops in media 
with initially very weak magnetic fields. The simulations follow the growth of magnetic 
fields and reproduce the characteristics of turbulence. These results are also compared with 
full cluster simulations that have examined the properties of ICM turbulence. 
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1. Introduction 

Observation and theory have revealed intra- 
cluster media (ICMs) to be very dynamic en- 
vironments with active "weather" driven by 
a host of activities such as mergers, accre- 
tion, AGNs, galactic winds and instabilities. 
These drivers are common and cause the ICMs 
to be criss-crossed by large-scaled, complex 
flows that generate shocks, contact discon- 
tinuities (aka "cold fronts") and bulk shear. 
Inevitably such flows should lead to turbulence 
in the ICMs, an outcome supported by grow- 
ing observational evidence. These include, for 
example, substantial ICM random velocities 
in Perseus reducing resonance scattering in 
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the 6.7 keV iron line IChurazov etal 2004 * ev- 
idence for thermal ICM pressure fluctuations 
in the Coma cluster lShuecker etal 20041 patchy 
Faraday rotation measure distributions in sev- 
eral clusters Bonafede etal 2010 and the ab- 
sence of large scale polarization in cluster ra- 
dio halos lKim etal 1990t suggesting disordered 
magnetic fields. 

Turbulence in clusters is important to un- 
derstand for many reasons. Turbulent pressure 
helps support the ICM, so relevant cluster mass 
measures. Turbulence transports entropy, met- 
als and cosmic rays, all important cluster di- 
agnostics. It transports and amplifies magnetic 
fields, which in turn control ICM viscosity, 
resistivity and thermal conductivity, as well 
as the propagation and acceleration of cosmic 
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rays. The literature on turbulence is extensive 
including excellent reviews on MHD turbu- 
lence, which is most relevant to the ICM (e.g., 
Brandenburg & Subramanian 2005| l. Here we 



make a few observations pertinent to this meet- 
ing. 

2. Origins of Cluster Turbulence 

Turbulence describes motions possessing sig- 
nificant random velocities. This random veloc- 
ity can include both compressible (V ■ « 0) 
and vortical, or solenoidal {a> = V x m 
0) components. In ICM circumstances, which 
usually involve mostly subsonic flows, the vor- 
tical component ordinarily dominates. Thus, 
understanding this turbulence begins with an 
identification of the sources of vorticity and the 
manner in which vorticity evolves. 

Euler's equation for an ideal fluid 
can be expressed in terms of vorticity 
ILandau & Lifshitz 19871 

^ = V X (m X -I- XrVp X 'VP. (1) 
Of 

This can be rewritten as a conservation law 
for the integrated vorticity through a surface 
element moving with the fluid, or by way of 
Stokes' Theorem a conservation of circulation 
around that surface element. 



applied to magnetic fields means that stretched 
magnetic flux tubes enhance local fields. The 
local B of an incompressible flux tube varies 
in proportion to the length of the tube. Flux 
tube stretching is, in fact, at the core of the tur- 
bulence dynamo, or small scale dynamo that 
amplifies weak magnetic fields inside turbu- 
lent conducting media. Vortical turbulent mo- 
tions statistically extend the length of a fluid 
element over time, causing both vorticity and 
magnetic field to be locally enhanced while 
conserving total circulation and magnetic flux 
in an ideal fluid. The root-mean-square (RMS) 
values of both vorticity and magnetic field in- 
tensity grow in this way, even in the absence of 
source terms. 

On the other hand, vorticity can be 
generated at the curved surface of shocks 
in and around clusters (e.g., |Ryu etal 2003"! 
Kang etal 20071 IVazza etal 2009l l and bv~gie 



baroclinity of flows. For uniform upstream 
flow, the vorticity behind a curved shock sur- 
face is 



(P2 - Pi)" j^rj 

Wcs = — KUi X n, 

PiPi 



(3) 



^ J co da = ^ ^u dl = J -^(VpxVP)-da, 

where d/dt is the Lagrangian time derivative. 

Equation |2] shows that in ideal, barotropic 
flows, where the pressure depends only on 
density (Vp X Vf = 0), a fluid element's 
net vorticity is conserved. Local values of ai 
still change, of course, without vorticity source 
terms. For instance, when an incompressible 
flow element's cross section is decreased as 
it is stretched, (o will increase. This vortex 
stretching (how tornados form) is central to 
evolution of turbulent flows. 

These vorticity properties are analogous to 
magnetic flux conservation familiar to all as- 
trophysicists. In fact, the magnetic induction 
equation in ideal MHD is the same as the 
barotropic form of equation [1] with the substi- 
tution cj -4 B. The vortex stretching analogy 



with pi and p2 the upstream and downstream 
gas densities, Ui the upstream flow velocity 
in the shock rest frame, K the shock surface 
curvature tensor, and h the shock normal unit 
vector If isopycnic (constant density) and iso- 
(2B)aric surfaces are not coincident, vorticity is 
generated according to equation |2] If we let 
P - Sp'^, where 5 is a proxy for entropy in 
a y-law gas, we see that the baroclinity is in- 
troduced by shock induced entropy variations. 
See |Ryu etal 2008 for more discussion of this 
issue. AGNs and galactic winds also add shear 
to ICMs, so equivalently, vorticity (w ~ Aujd, 
where 6 is the thickness of the shear layer). 

Subsequently, the energy in the vortical 
flows cascades down to smaller scales and tur- 
bulence develops, provided the viscous dissi- 
pation scale, U.isc, is smaller than the "driving 
scale", Ld, for the flow. The driving scale, L,/, 
is generally comparable to such things as the 
curvature radius of a shock, the size of a sub- 
structure core, or the scale of an AGN or galac- 
tic outflow. These likely range from ~ 10s of 
kpc to ~ 1 00s of kpc. 
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The appropriate viscous dissipation scales 
in ICMs remain uncertain. The media are hot 
and very diffuse plasmas, so Coulomb colli- 
sions are ineffective. The associated mean free 
path, Acoui ~ 1 kpc T^^v/(n_3Uth,ioo), ranges 
from ~ 10s of pc to ~ 10s of kpc, depend- 
ing on the cluster circumstances. Here Tkgv is 
the ICM temperature in keV, n_3 is the density 
in lO^^cm"^ and m,/,j()o is the thermal veloc- 
ity of ions in 100 km/sec. The corresponding 
kinetic viscosity, v ~ Uih'^Coui, is very large, 
and the Reynolds number. Re ~ LdU/v ~ 
few X 10 , with U ~ u,/, ~ few to 10s, the 
velocity at the driving scale. Hence, the vis- 
cous dissipation scale due to Coulomb scatter- 
ing alone, ~ Ld/Re^^*, would range from 
fractions of a kpc in cooling cores to several 
10s of kpc in hot ICMs. On the other hand, 
it has been suggested that plasmas threaded 
with weak magnetic fields, like the ones in the 
ICMs, are subject to gyro-scale instabilities, 
such firehose and mirror instabilities; then, 
the scattering of particles with the resulting 
fluctuations could reduce the particle mean- 
free path, so also the viscous dissipation scale 



( [Schekochihin & Cowley 2006| l. The detailed 
picture is still very uncertain. Nevertheless, in 
our discussion below we assume the physical 
dissipation scale is at least as small as the ef- 
fective dissipation scales of the simulations, so 
of order the grid resolution. In our simulations, 
the grid resolution would be ~ 1 kpc or so (see 
§3), while for full cluster simulations it is at 
least several kpc. 

The resistive dissipation scales. Ires, in the 
ICMs are also uncertain. In a turbulent flow 
with ?7 «; V, it is Ires ~ Ld/Rm^^^, where Rm ~ 
LjU/r] is the magnetic Reynolds number, and 
J] is the resistivity. It is likely that the magnetic 
Prandtl number P, ^ = Rni/Re = v/77 ^ 1 
in the ICMs. For instance, viscosity and re- 
sistivity due exclusively to Coulomb scattering 
would lead to P^.m » 1 (e.g., Spitz er 1978| . 
In the simulations reported here, both the vis- 
cosity and resistivity have numerical origins 
with dissipation scales of order the grid reso- 
lution; thus, Pr^M ~ 1- Most simulations, in- 
cluding most full cluster simulations, effec- 
tively have Pr^M ~ 1 too. Strictly speaking, 
the simulations with f r,M ~ 1 should be valid 



only in, and so applied only to the scales of 
/ ^ and I ^ For the scales of Ires ^ 
/ ^ /,.,jc, simulations with large Prandtl num- 
ber are required. It is, however, hard to re- 
produce large Prandtl number turbulence with 
sufficient viscous and resistive inertial ranges 
(e.g. . .Schekochihin etal 2004J . 

3. Simulation of ICI\1-like Turbulence 

It is not simple to isolate turbulence from gen- 
erally complex motions in real or simulated 
clusters. As a complementary exercise we have 
initiated a high-resolution simulation study of 
the evolution and saturation of driven MHD 
turbulence in periodic computational domains 
that resemble ICMs. Since cluster media, while 
clearly magnetized, are not energetically dom- 
inated by magnetic fields, we focus on turbu- 
lence developed with initially very weak mag- 
netic fields. The full study considers both com- 
pressible and incompressible fluids and ideal 
and nonideal media with a range of magnetic 
Prandtl numbers. We report here some ini- 
tial results of simulations of ideal, compress- 
ible MHD turbulence in isothermal media. The 
simulations used an isothermal ideal MHD 
code, which is a version updated from that 
of IKim etal 20011 for massive parallelization. 
Results here are from simulations with 1024-' 
and 2048^ grid zones (with the typical cluster 
sizes ~ 1 Mpc, the grid resolution would be 
~ 1 kpc). The two simulations are very consis- 
tent; the higher resolution calculation produces 
slightly stronger magnetic fields at saturation. 

Initially the medium has uniform density, 
p = 1, with gas pressure, P - 1, so that 
the sound speed, c, = 1. The initially mag- 
netic field is uniform and very weak, with /3 - 
P/Pb = 10''. The cubic box has dimensions, 
Lo = = Ly = L- = 10 with periodic bound- 
aries. Turbulence is driven by solenoidal ve- 
locity forcing, drawn from a Gaussian random 
field determined with a power spectrum P^ oc 
k^exp(-Sk/kpeak), where kpeak = 2^o (^0 = 
IjtILq), and added at every Af = O.OIL/Cj. 
The power spectrum peaks around ~ I .Sko, 
or around a scale, Ld 2/3Lo- The am- 
plitude of the perturbations is tuned so that 
urms ~ 0.5 or Ms = urms/cs ~ 0.5 at satu- 
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ration, close to what resuilts in full cluster sim- 
ulations (e.g., |Nagai etal 2007| [Ryu"etal 2008 



IVazza etal 20091 and see also §4). 

Our setup gives a characteristic time scale, 
td — LdjuRMS ~ 15. In that time the largest 
eddies should spawn something resembling 
turbulent behavior. Indeed Fig. [T] shows that 
the mean turbulent kinetic energy density, Ek, 
grows and peaks at f ~ f^, with a value cor- 
responding to urms ~ 2/3. Subsequently, Ek 
slowly declines as the mean magnetic energy 
density, Eb grows. The kinetic energy power 
spectrum, E^ik), at f = 20 also shown in Fig. 
[1] exhibits a peak at k/ko ~ 2, near the driv- 
ing scale, and takes a Kolmogorov-like, iner- 
tial form, E(k) oc k-^/\ for k/ko ^ 50. By 
this time energy has cascaded from the driving 
scales far enough that the motions, with neg- 
ligible magnetic backreaction, are reasonably 
described as classical, hydrodynamical turbu- 
lence over a modest range of scales. 

Turbulent fluid motions stretch and inten- 
sify vorticity and magnetic flux, leading to de- 
velopment of a chaotic sea of vortex and mag- 
netic filaments. This is illustrated for the mag- 
netic field at f = 20 in the top of Fig. |2] The 
magnetic filaments in this image have charac- 
teristic lengths of a few % of Lq. This agrees 
with the fact that the magnetic power spectrum, 
Esik), peaks for kpeaklko ~ 50 at this time. The 
transverse dimensions of the filaments, an or- 
der of magnitude smaller at this time, seem to 
correspond to the dissipation scale. 

As the magnetic field intensifies, both vor- 
ticity and magnetic energies undergo inverse 
cascades from small to large scales, with the 
coherence lengths of their filaments growing 
accordingly. This is evident for the magnetic 
field in the power spectrum changes in Fig. 
[T] The inverse cascade of magnetic energy 
can be understood as follows. The field is 
wrapped more quickly around smaller scale ed- 
dies, because the eddy turn over time varies 
as f/ oc f^^. Maxwell stresses, oc (V x Z?) x 
B, then, feed back on the kinetic turbulence, 
causing significant modifications in the fluid 
motions, thus saturating the magnetic field 
growth on a given scale, /, when £^(0 ~ 
Efcil)- Since the turbulent kinetic energy on 
a scale Ek{1) oc f-^^^ the saturation scale of 
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Fig. 1. Top: Evolution of kinetic, Ek, and mag- 
netic, Eb, energies in simulations of ideal, 
compressible MHD turbulence for two grid 
resolutions. Bottom: Power spectra, E{k), of 
kinetic and magnetic energies at f = 20 and 
t — 130 in the 2048^ zone turbulence. 



the magnetic turbulence should evolve over 
time as Ib <^ t"^^, while the magnetic energy 
grows linearly with time, both consistent with 
Fig.[T] Eventually, as Ib approaches the outer 
scale of the kinetic turbulence, Ld, the scal- 
ings break down and turbulence reaches sat- 
uration where the ratio of the total magnetic 



Jones: Cluster Tubulence 



5 



'i- 




Fig. 2. Magnetic energy density distributions in 
MHD turbulence. Top: Log(£'B) at f = 20 in 
2048^ simulation. Bottom: Log(£'B) at t- 130. 
Higher tones at the late time reflect stronger 
fields. 



to kinetic energy is EbIEk ~ ll'i (see also, 
e.g., ISchekochihin etal 20041 |Ryuetal2008 



ICho etal 2009"! |Cho & Ryu 2009f 



Finally, we emphasize an interesting topo- 
logical transformation in the flow structure as 
turbulence proceeds through the linear growth 
to the saturation stage. Fig. |2] displays the dif- 
ferent topologies of the magnetic flux struc- 



tures at f = 20 and f = 130. At the ear- 
lier time the field is organized into individ- 
ual filaments. At the later time those filaments 
have evolved into striated, ribbon-like forms 
(see also ISchekochihin etal 2004l l. Close ex- 
amination reveals the ribbons to be laminated, 
with magnetic field and vorticity interleaved 
through each cross section on scales of the or- 
der the dissipation length. In hydrodynamical 
turbulence such ribbons would be unstable, but 
the interplay of vorticity and magnetic field 
seems to stabilize them in MHD turbulence. 

The t - 130 image in Fig. |2] also high- 
lights the important fact that the magnetic 
field in MHD tubulence is highly intermittent. 
Relatively strong field ribbons wrap around 
large shear leaving weak field cavities. Such 
structures are quite distinct from what one ob- 
tains, for example, if they construct a mag- 
netic field distribution out of a Gaussian ran- 
dom variate, even if the outcome is a mag- 
netic field distribution having exactly the same 
power spectrum, as illustrated very nicely in 
previous work bv iWaelkens etal 20091 

4. Comparison to Cluster Simulations 

Recent high-resolution cluster formation sim- 
ulations have been analyzed for informa- 
tion about properties of ICM turbulence 
and its evolution. A few of these simula- 
tions include MHD e.g., IDonnert etal 20091 
IRuszkowski etal 20101 IXu etal 20101 With or 
without magnetic fields, an initial challenge 
is identifying true turbulence in generally 
complex flows, especially during and after 
merging activity. Purely solenoidal motions 
in simulations in a periodic domain can be 
cleanly isolated using Fourier transforms (e.g., 
Ryu etal 2008| l. Otherwise some kind of spa- 



tial filtering is needed that analyses the motions 
only up to a maximum scale, such as the core 
radius of the cluster (e.g., |Doiag etal 2005 ). We 
mention here a few findings of general signifi- 
cance in this context. 

Several studies have emphasized turbu- 
lence generation from shocks in and around 
clusters, especially coming from structure for- 
mation generally, and merger activity specif- 
ically (e.g., |Kang etal 2007] |Ryu etal 2008 



6 



Jones: Cluster Tubulence 



IVazza etal 2009l l. Several studies found that 
turbulent energy in post-merger ICMs can 
be comparable to, although somewhat smaller 
than, the thermal energy; it commonly reaches 
levels Ek ~ 1 / 10 - 1/4 Ej (|Nagai etal 2007 
Ryu etal 2008] IVazza etal 20091 1, similar to 



our simulations discussed in §3. Thus, tur- 
bulent pressure may impact on hydrostatic 
equilibrium-based mass estimates. Simulations 
also suggest that the relative turbulent pres- 
sure support is greater towards cluster out- 
skirts (e.g.,|Ryu etal 2008||Lauetal 200911. The 



turbulence formed in full cluster simulations 
does not have a sufficiently wide inertial range 
to evolve a true Kolmogorov power spec- 
trum. Within that limitation, however, the 
reported kinetic energy power spectra are 
consistent with expectations. Several studies 
demonstrated that cluster turbulence can am- 
plify magnetic fields to at least /zGauss levels 
(e.g., IDonnert etal 20091 IXu etal 20ld] l. This 
amounts to ~a % or so of the thermal pres- 
sure and ~10 % or so of the kinetic turbu- 
lent pressure. From simulations such as we 
reported here, the time to reach equipartition 
(Eb ~ 2/3Ek) from an initially weak field is 
many driving scale times. It is not surprising 
that the fields produced in clusters are well be- 
low those levels. By the same token the mag- 
netic field power spectra may peak well below 
the driving scales. The magnetic power spec- 
trum reported from cluster formation simula- 
tions using magnetic fields seeded by AGNs in 
IXu etal 20101 for example, is consistent with 
this expectation. 

Turbulence decay generally takes several 
eddy times on the driving scale once the driv- 
ing ends. In clusters we expect decay times, 
Td ^ L^/.ioo/AwiooGyr, where Ljjoo in 100 kpc. 
This is consistent with turbulence decay times 
~ Gyr reported in cluster formation simula- 
tions (e.g.. lPaul etal 2010l l 

5. Conclusions 

Processes such as shocks and outflows are 
likely to drive turbulence in ICMs. The detailed 
physics is difficult to model analytically, but 
simulations allow us to explore it in some de- 
tail. Magnetic fields are integral ingredients of 



both the microphysics of ICM transport prop- 
erties and essential players in the large scale 
dynamics. Simulations are revealing important 
insights into the character of the turbulence, in- 
cluding properties of magnetic fields. 
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